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Abstract. The coherent and collective nature of Bose-Einstein condensate can en- 
hance or suppress physical processes. Bosonic stimulation enhances scattering in al- 
ready occupied states which leads to atom amplification, and the suppression of dis- 
sipation leads to superfluidity. In this paper, we review several experiments where 
£j ' suppression and enhancement have been observed and discuss the common roots of 

and differences between these phenomena. 
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^ ■ When a gas of bosonic atoms is cooled below the transition temperature of Bose- 

Einstein condensation, it profoundly changes its properties. The appearance of a 
macroscopically occupied quantum state leads to a variety of new phenomena which 
set quantum fluids apart from all other substances. Fritz London even called them 
the fourth state of matter [1]. 

Many of the key concepts in quantum fluids were derived from studying the 
weakly interacting Bose gas, for which rigorous theoretical treatments were possible 
[2,3]. In 1995, with the discovery of BEC in a dilute gas of alkali atoms [4-6], it 
became possible to study such a system experimentally . The theoretical framework 
connects the observed equilibrium and dynamic properties to the presence of long- 
range order and low- lying collective excitations [7]. 

Many special properties of Bose condensates involve the suppression or enhance- 
ment of physical processes. Our recent experiments include the suppression and 
enhancement of elastic collisions of impurity atoms [8], the suppression of dissi- 
pation due to superfluidity [9,10], and the suppression [11] and enhancement of 
light scattering [12]. Bosonically enhanced Rayleigh scattering was used to amplify 
either atoms [13] or light [14] in a condensate dressed by laser light. We review 
these experiments and discuss the properties of the Bose condensate which lead to 
enhancement and suppression. 
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I SCATTERING OF LIGHT AND MASSIVE 

PARTICLES 



Before we discuss light scattering and collisions in a BEC, we want to derive 
some simple general expressions based on Fermi's golden rule which will be useful 
to see the similarities and differences between the different processes. 

When a condensate scatters a photon or material particle, the scattering is de- 
scribed by the Hamiltonian 

H' — C c\a) n c k a m 5 l+n _ k _ m (f) 

k,l,m,n 

Here c k (c k ) is the destruction (creation) operator for the scattered particles, and 
a k (a k ) is the destruction (creation) operator for atomic plane waves of wavevector 
k. The strength of the coupling is parametrized by the coefficient C. 

We consider the scattering process where a system with N atoms in the conden- 
sate ground state and N q quasi-particles with wavevector q scatters particles with 
incident momentum k into a state with momentum k — q. The initial and final 
states are 

\i) = \n k ,n k _ q ;N ,N q } 

\f) = \n k - 1, n k ^ q + 1; N - 1, N q + 1) (2) 

respectively, where n k denotes the population of scattering particles with wavevec- 
tor k 1 . It should be emphasized that, due to the interatomic interactions, the 
quasi-particles with occupation N q are not the plane waves created by the operator 
a q , but the quanta of collective excitations with wavevector q. 

The square of the matrix element Mi between the initial and final state is 

IM^HUW)! 2 

= \C\ 2 \(N — N — 1, N q — l\fi(q)\N = N,N q = 0}\ 2 (N q + \){n k _ q + l)n k (3) 

where p(q) = J2 m aL+^m is the Fourier transform of the atomic density operator 
at wavevector q. 

The static structure factor of the condensate is S(q) = (g\p{q.) ft {<l)\g) / N where 
\g) = \N = N, N q = 0) is the BEC ground state. We then obtain for the scattering 
matrix element Mi 

(Mil 2 = \C\ 2 S(q)(N q + l)(n fc _, + 1 )N n k (4) 

The scattering rate W\ for the process \n k , n k _ q ; N , N q ) — > \n k — l,n k _ q + l; N — 
1, N q + f ) follows from Fermi's golden rule as 

^ This choice of final states implies that we neglect scattering between quasi-particles and con- 
sider only processes involving the macroscopically occupied zero-momentum state of the conden- 
sate. Formally, we replace the Hamiltonian (Eq. 1) by C^2 k q {cl_ q a q dka + c* k _ q a\cka- q ). 




the energy of quasi-particles with momentum q (which we will later obtain from 
Bogoliubov theory). To obtain the net scattering rate, one has to include the reverse 
process \n k ,n k _ q ; N , N q ) — > \n k + l,n k _ q — l; N + l, N q — 1) by which atoms scatter 
back into the condensate. The square of the matrix element M 2 for this process is 
\C\ 2 S(q)N q n k - q (No + l)(rik + 1). The net rate W + of scattering atoms from the 
condensate into the quasi-particle mode q is the difference of the two partial rates 
W + = Wi — W 2 . Assuming N 3> 1 (i.e. N + 1 m N Q ), we obtain for the net rate 



W + = — \C\ 2 S(q)N [n k (N q + n k _ q + 1) - N q n k . q ]6(E k - E k _ q - ftu*) (6) 



stimulation term (N q + n k _ q + l) is approximately (JV + 1) with N = max(N q , n k - q ). 
This illustrates that there is no bosonic stimulation of the net rate by the least 
populated final state. With the dynamic structure factor S(q,u>) = S(q)S(u — uj^) 
Eq. 6 simplifies to 



The rate W+ in Eq. 6 is the rate for the Stokes process where E k > E k _ q . 
Momentum transfer q to the condensate is also possible as an anti-Stokes process 
where a quasi-particle with momentum — q is scattered into the condensate, and 
the scattered particle gains energy. The net rate W_ for this process is obtained in 
an analogous way as 



W. = —\C\ 2 S(q)N [n k (N_ q - n k „ q ) - (N. q + l)n k _ q }5{E k - E k . q + huf). (8) 



The net scattering rates in Eqs. 6 and 8 are the product of three terms. 

• The static structure factor S(q) represents the squared matrix element for the 
condensate to absorb momentum q. 

• The delta function denotes the density of final states. 

• The bosonic stimulation term represents stimulation by the occupation in the 
final state either of the scattering particles or the condensate. 

The interplay of these three terms is responsible for the suppression and enhance- 
ment of physical processes in a condensate. The properties of the condensate as 
an intriguing many-body system are reflected in the structure factor. In section 
II, we discuss its measurement through stimulated light scattering. The density 
of states is responsible for superfluidity because it vanishes for initial velocities of 
the incident particles which are smaller than the Landau critical velocity (Sections 
III and IV). Finally, bosonic stimulation by the occupancy N q of final states was 
responsible for superradiance, matter wave amplification and optical amplification 
in a condensate (Section V). 
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FIGURE 1. Bragg scattering of phonons and of free particles. Momentum transfer per particle, 
in units of fiq, is shown vs. the frequency difference uj/2-k between the two Bragg beams. The 
two spectra were taken at different atomic densities. The open symbols represent the phonon 
excitation spectrum for a trapped condensate (at a chemical potential fi/h = 9.2 kHz, much 
larger than the free recoil shift of ~ 1.4 kHz). Closed symbols show the free-particle response of 
a twenty-three times more dilute (ballistically expanded) cloud. Lines are fits to the difference 
of two Gaussian line shapes representing excitation in the forward and backward directions. See 
Refs. [11,16] for more details. Figure is taken from Ref. [11]. 

II DETERMINATION OF THE STATIC STRUCTURE 
FACTOR BY BRAGG SPECTROSCOPY 

The BEC matrix element or the structure factor in Eqs. 6 and 8 can be directly 
determined experimentally by stimulated light scattering. The density of states 
(the 5 function in Eq. 6) does not restrict the scattering since the photon energy is 
much higher than the quasi-particle energy. As a result, photons can be scattered 
into the full solid angle and provide the necessary recoil energy of the atom by a 
small change (~ 10~ 9 ) in the frequency of the scattered photon. 

When the condensate is illuminated by two strong laser beams with wavevectors 
k and k — q and a difference frequency u, the rate W of transferring photons from 
one beam to the other is 

W/N = (W+ + W-)/N = 2iruj R S(q)(5(uj - cuf) - 5(u + uf)) (9) 

where the two-photon Rabi frequency uj r is given by uj 2 r = \C\ 2 rikn k _ q /Ti 2 . When 
D is scanned the "spectrum" of a condensate consists of two peaks at ±o>,f (Fig. 
1). The strength (or integral) of each peak corresponds to S(q). We refer to this 
method as Bragg spectroscopy since the basic process is Bragg scattering of atoms 
from an optical standing wave [15]. 

We observed that the scattering rate was strongly suppressed when the momen- 
tum transfer frq became smaller than the speed of sound c in the condensate, i.e., 
when the light scattering excited a phonon and not a free particle. These observa- 
tions are in agreement with the Bogoliubov theory, which obtains the elementary 
excitations or Bogoliubov quasiparticles as eigenstates of the Hamiltonian (Eq. 
1). A quasi-particle with wavevector q is annihilated by the Bogoliubov operator 



b q = u q a q + v q CL- q , where u q = cosh0 g , v q = sinh0 9 and tanh20 g = fj,/(hu) q + fi). 
Its energy is given by 

hu* = v / h o (H° + 2 ^)' ( 10 ) 

where hu q = h 2 q 2 /2m and fi is condensate's chemical potential. Expressing p(q) by 
Bogoliubov operators leads to an expression for the static structure factor S(q) = 
(u q — Vq) 2 . The static structure factor tends to S(q) — > hq/2mc and vanishes in 
the long wavelength limit, as required of a zero-temperature system with finite 
compressibility [17]. 

This reveals a remarkable many-body effect. For free particles, the matrix ele- 
ment for momentum transfer is always 1, which reflects the fact that the operator 
e iq ' r connects an initial state with momentum p to a state with momentum p + hc\ 
with unity overlap. For an interacting Bose-Einstein condensate, this overlap van- 
ishes in the long- wavelength limit. As we have discussed in previous publications 
[11,16], this suppression of momentum transfer is due to an destructive interference 
between the two pathways for a condensate to absorb momentum q and create 
a quasi-particle: one pathway annihilates an admixture with momentum — q, the 
other creates an admixture at momentum +q. 

The leading terms in stimulated light scattering do not depend on temperature. 
The rates W + and W- (Eqs. 6 and 8) are independent of the thermally excited 
population of quasiparticles N q and N_ q in the limit of large n k , n k - q ^> N± q , i.e. 
when the scattering is induced by a second laser beam. For spontaneous scattering 
(n k ^ q = 0), one obtains for the total scattering rate W instead of Eq. 9 

2lT r i 

W/N = ^\C\ 2 S(q)n k [(1 + N q )(6(u - uf) + N_ g 5(cu + uf))] . (11) 

Absorption of and bosonic stimulation by thermally excited quasi-particles become 
important when the temperature is comparable to or larger than the quasi-particle 
energy hujq. Eq. 11 is proportional to the temperature dependent dynamic struc- 
ture factor Sr(q, w) [18]. 



Ill SUPPRESSION OF IMPURITY COLLISIONS 

The key difference between the scattering of light and massive particles (or im- 
purities) is their energy-momentum dispersion relation. The dispersion relation for 
impurities is E k = (hk) 2 /2M, whereas for light, E k = hkQ with q denoting the 
speed of light. This difference is responsible for the complete suppression of impu- 
rity scattering at low velocities. For simplicity, we assume equal mass M for the 
impurity and condensate atoms. 

Energy-momentum conservation, the delta function in Eq. 6, requires E k —E k _ q = 
hujq . For impurity particles, the l.h.s. is always less than vhq, where v = hk/M is 
the initial velocity of the impurities. Thus, collisions with the condensate are only 
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Impurity velocity [in units of the speed of sound] 
FIGURE 2. Suppression of collisions. Shown is the suppression factor as a function of the 
impurity velocity (normalized by the speed of sound c, solid line). The dotted line represents the 
suppression due to phase-space restriction alone (i.e., setting the structure factor S(q) = 1). 



possible, if this maximum energy transfer is sufficient to excite a quasi-particle, i.e., 
v > min(ujq/q) = Vl, where vl is the Landau critical velocity [19] for superfluidity 
below which no dissipation occurs because the density of final states vanishes. 
For the excitation spectrum of the condensate (Eq. 10), the Landau velocity 

is the Bogoliubov speed of sound vl = c = yJ/i/M. Impurity particles moving 
below this speed cannot dissipate energy in collisions. In contrast, for photons, 
dEk/d{hk) = a ^> c, and scattering is always possible. 

In the perturbative limit (no stimulation by the final occupation), the total rate 
of scattering T is given by integrating Eq. 6 over all possible momentum transfers. 
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= (JVb/V) a(rj) v, (12) 

where HQ = Mv(l — l/i] 2 ) is the maximum possible momentum transfer, and 
i] = v/c must be larger than 1. We have used C = Aixh 2 a/MV for s-wave scatter- 
ing between condensate and the impurity particles where a denotes the scattering 
length and V the condensate volume. The collision cross section is cr(rj) = aoF(rj) 
where <7n = 47ra 2 . For r) < 1, F(rf) = and for rj > 1, F(rj) = 1 — l/r/ 4 — log(r] 4 )/?7 2 . 

The suppression factor ^(77) is determined by two factors: the phase space re- 
striction due to the delta function in Eq. 12, and additional suppression at low 
momentum transfers by the structure factor of the condensate. For decreasing ve- 
locity, the possible scattering angles become restricted to a forward scattering cone, 
which shrinks to zero solid angle at the Landau critical velocity. A graph of the 
suppression factor as function of impurity velocity is shown in Fig. 2. 
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FIGURE 3. Observation of elastic collisions between the condensate and impurity atoms, (a) 
Impurities traveling at 6 cm/s along the radial axis (to the left in images) were scattered into a 
s-wave halo. Absorption image after 50 ms of time-of-flight shows the velocity distribution after 
collisions between the condensate (bottom) and the outcoupled tuf = atoms (top), spatially 
separated by a Stern-Gerlach type magnetic field gradient. The collisional products are distributed 
over a sphere in momentum space. The image is 4.5 x 7.2 mm. (b) Similar image as (a) shows the 
collisional products (arrow) for impurity atoms (top) traveling at 7 mm/s along the condensate 
axis (upward in image). For this image, v g /c = 2.7 (see text). Collisions are visible below the 
unscattered impurities, (c) Similar image as (b) with v g /c = 1.6. Collisions are suppressed. The 
outcoupled atoms were distorted by mean-field repulsion. The images are 2.0 x 4.0 mm. 

Experimentally, the Landau critical velocity can usually only be observed by 
moving microscopic particles through the superfluid which do not create a macro- 
scopic flow pattern. Studies of superfluidity with microscopic objects were pursued 
in liquid 4 He by dragging negative ions through pressurized 4 He [20,21], and by 
scattering 3 He atoms off superfluid 4 He droplets [22]. 

To study the effects of impurities interacting with the condensate, we created 
microscopic impurity atoms using a stimulated Raman process which transferred 
a small fraction of the condensate atoms in the \F = l,m^ = —1) hyperfine state 
into an untrapped hyperfine state \F = l,m,F = 0) with a well-defined initial 
velocity [23]. The initial velocity could be adjusted between zero and two recoil 
velocities by varying the angle between the two Raman beams. As these impurities 
traversed the condensate, they collided with the stationary condensate, resulting 
in a redistribution of the impurity momenta which was detected by a time-of-flight 
analysis [8] (see Fig. 3). 

To probe for the suppression of collisions, one has to vary the impurity velocity 
around the speed of sound. For that, we produced impurity atoms at low velocities 
(7 mm/s) and varied the speed of sound by changing the condensate density. The 
small axial velocity imparted by Raman scattering allowed us to identify products 
of elastic collisions in time-of-flight images (Fig. 3b, c) since collisions with the sta- 
tionary condensate redistributed the impurity atoms toward lower axial velocities. 
However, the impurity velocity was predominantly determined by the gravitational 
acceleration g, which imparted an average velocity of v g = \/2gz c where z c is the 
Thomas- Fermi radius of the condensate in the z-direction. Thus, the effect of su- 
perfluidity on impurity scattering depends primarily on the parameter fj = v g /c 
which is the ratio of the typical impurity velocity v g to the speed of sound c at the 
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FIGURE 4. Onset of superfluid suppression of collisions. Shown is the observed collisional 
density normalized to the predicted one in the limit of high velocities as a function of rj = v g /c, 
which is a measure of the impurity velocity in units of the condensate's speed of sound. The error 
bars represent the statistical uncertainty. Data is taken from Ref. [8] 
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FIGURE 5. Collective amplified elastic scattering in a Bose-Einstein condensate. Shown is the 
fraction of collided atoms vs. the number of outcoupled atoms. For this data, w g /c=4.9 and the 
chemical potential was 1.8 kHz. Figure is taken from Ref. [8] 



center of the condensate. 

A time-of-flight analysis of impurity scattering for the case of a low-density con- 
densate (small c) and large condensate radius (large v g ) is shown in Fig. 3b. The 
effect of collisions is clearly visible with about 20% of the atoms scattered to lower 
axial velocities (below the unscattered impurities in the image). In contrast, in the 
case of tight confinement, the condensate density is higher (larger c) and its radius 
is smaller (smaller v g ), and the collision probability is greatly suppressed due to 
superfluidity (Figs. 3c and 4). 

The bosonic stimulation factor in Eq. 6 is relevant if the final states are populated, 
either by scattering or thermally. We observed that the fraction of collided atoms 
increased with the number of outcoupled impurities (Fig. 5). For a large outcoupled 
fraction, population n^_ q and N q is built up in the final states and stimulates further 
scattering. This collisional amplification is not directional, and is similar to the 
recently observed optical omnidirectional superfluoresence [24]. 

Gain of momentum and thus transfer of energy from the condensate to the impu- 
rity atoms is impossible at zero temperature, but may happen at finite temperature 
due to the presence of thermal excitations (the N_ q term in Eq. 8) 2 . Thus finite 



2 > We are grateful to S. Stringari for pointing out the importance of finite-temperature effects. 
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FIGURE 6. Temperature dependent cross-section vs. impurity velocity. Shown is the 
cross-section at zero temperature (solid line) and at a finite temperature kT = fj, which is typical 
for our experimental conditions (dash-dotted line) . The finite temperature cross-section includes 
collisions involving thermally occupied quasi-particles where the impurities lose or gain energy. 
In the experiment, we measured the number of impurities which lost its energy minus the number 
which gained energy. Thus, the experimental measured cross-sections (Fig. 4) should be compared 
to <7 co ii,ioss - fcoii.gain (dashed line). 

temperature enhances the elastic cross section by two effects: Absorption of quasi- 
particles (anti-Stokes process) and stimulation of momentum transfer by the final 
state population (Stokes process). 

Fig. 6 shows the dramatic variation of the elastic scattering cross section with 
temperature. However, the finite temperature did not affect our data in a ma- 
jor way: Due to gravitational acceleration we couldn't probe the velocity regime 
well below the Landau critical velocity where only thermally assisted collisions are 
possible. Furthermore, when we counted the number of collided atoms we had to 
use a background subtraction method where we subtracted the small signal of the 
energy gain collisions from the energy loss collisions (see Ref. [8] for details), thus 
cancelling most of the finite-temperature effects. 

IV SUPPRESSION OF DISSIPATION FOR A MOVING 

MACROSCOPIC OBJECT 

So far, we have discussed the suppression and enhancement of microscopic pro- 
cesses (light scattering and impurity collisions). The suppression of dissipation is 
even more dramatic on the macroscopic scale. The flow of liquid 4 He and the mo- 
tion of macroscopic objects through it are frictionless below a critical velocity [25]. 
Recently, we have explored such frictionless flow in a gaseous BEC [9,10]. 

The microscopic and macroscopic cases bear many parallels. The onset of scat- 
tering or dissipation has two requirements: one needs final states which conserve en- 
ergy and momentum, and an overlap matrix element which populates these states. 
In the case of macroscopic flow, the first requirement leads to a critical velocity 



for vortex creation and the second requirement addresses the nucleation process of 
vortices. 

The Landau criterion for superfluidity shows that excitations with momentum p 
and energy E(p) are only possible when the relative velocity between the fluid and 
the walls or a macroscopic object exceeds the Landau critical velocity vl which is 
given by v L = mm(E(p)/p) (see e.g. [2,25]. A similar criterion applied to vortex 
formation yields 

-^vortex h ( d\ , s 

"^W-M^UJ <13) 

where I = J Pd 3 r is the integrated momentum of the vortex ring or line pair, E 
is its total energy, d is the dimension of the container, and £ the core radius of a 
vortex which in the case of dilute gases is the healing length £ = 1/ \f%ixna. Ref . [25] 
derived Eq. 13 for vortex rings with a maximum radius d, Ref. [26] looked at pairs 
of line vortices at distance d. Feynman [27] found a similar result for superflow 
through a channel of diameter d. 

A similar result is obtained for a Bose condensed system placed under uniform 
rotation with angular velocity Q. A vortex becomes energetically allowed when its 
energy E' in the rotating frame drops to zero, 

E' = E — VLL = (14) 

where E and L are the energy and angular momentum in the laboratory frame. 
This defines a critical angular velocity below which a vortex cannot be sustained 
due to conservation of angular momentum and energy 



The critical velocity at the wall of the rotating container, v c = dQ c , agrees with 
Eq. 13. 

However, Eqs. 13 and 15 only reflect the energy and momentum required to 
generate vortices, and do not take into account the nucleation process. If the 
scattering particle is macroscopic in size, the coupling is between the ground state 
and a state containing a vortex. Populating such a state requires nucleation of 
the vortex by the perturbing potential, which usually does not occur until higher 
velocities are reached than those predicted by Eqs. 13 and 15. The other option, the 
formation of the vortex by macroscopic quantum tunneling between the two states 
is an extremely slow process. In recent experiments in which a Bose condensate 
was placed in a rotating potential, the critical angular velocity for the formation 
of a single vortex was observed to be 1.6 times higher than the value given by Eq. 
15 [29]. This discrepancy may be due to a nucleation barrier associated with the 
excitation of surface modes, as some authors have recently suggested [30,31]. 



To study frictionless flow in a Bose-condensate, we focused a blue-detuned 514 
nm Argon laser beam onto the sodium condensate, which repelled atoms from the 
focus. The laser beam was scanned back and forth along the axial direction of the 
condensate, creating a moving "hole" that simulated a macroscopic object. Rapid 
sequence phase-contrast imaging allowed us to directly measure the flow pattern of 
the superfluid around the moving laser beam. 

For a weakly interacting Bose-condensed gas at density n(r) and chemical po- 
tential /i, pressure is identical to the mean- field energy density P = /i(r)n(r)/2 [7]. 
A drag force arises due to the pressure difference across the moving object. The 
chemical potential is given by fi(r, t) = gn(r, t), where g = Anh 2 a/M is the strength 
of two-body interactions. The drag force F is given by 

F ~ gSnAn = SfiAfi/g (16) 

where An and Afi are the differences in density and chemical potential across 
the stirring object, and S the surface area the macroscopic object presents to the 
condensate. 




FIGURE 7. Pressure difference across a laser beam moving through a condensate. On the 
left side in situ phase contrast images of the condensate are shown, strobed at each stirring half 
period: beam at rest (top); beam moving to the left (middle) and to the right (bottom). The 
profiles on the right are horizontal cuts through the center of the images. The stirring velocity 
and the maximum sound velocity were 3.0 mm/s and 6.5 mm/s respectively. Figure is taken from 
Rcf. [10]. 



If the laser beam is stationary, or moves slowly enough to preserve the superfluid 
state of the condensate, there will be no gradient in the chemical potential across 
the laser focus, and therefore zero force according to Eq. 16. A drag force between 
the moving beam and the condensate is indicated by an instantaneous density dis- 
tribution n(r, t) that is distorted asymmetrically with respect to the laser beam. 
Fig. 7 shows phase contrast images strobed at half the stirring period (where the 
laser beam is in the center of the condensate). A bow wave and stern wave form 
in front of and behind the moving laser beam, respectively. We define the asym- 
metry A as the relative difference between the peak column densities in front (72/) 



and behind (n h ) the laser beam A = 2(77/ — h b )/(h f + h b ). The asymmetry A is 
proportional to the drag force F. 

In Fig. 8 we show measurements of the asymmetry for two maximum densities 
n of 9 x 10 13 and 1.9 x 10 14 cm -3 . In each data set there is a threshold velocity 
v c below which the drag force is negligible, and this threshold increases at higher 
density. Its value is close to 0.1 c for both data sets, where c is the sound velocity. 
Above this critical velocity, the drag force increases monotonically, with a larger 
slope at low density. In addition, the heating rate due to friction against the laser 
beam was directly measured through time of flight absorption imaging [9,10,32] 
and found to be in good agreement with the value F • v, with F estimated from 
Eq. 16. 
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FIGURE 8. Density dependence of the critical velocity. The onset of the drag force is shown 
for two different condensate densities, corresponding to maximum sound velocities of 4.8 mm/s 
(•, left axis) and 7.0 mm/s (x, right axis). The stirring amplitudes are 29 ^m and 58 /im, 
respectively. The two vertical axes are offset for clarity. The bars represent statistical errors. 
Figure is taken from Ref. [10]. 



The observed critical velocity may be related to the formation of vortices. An 
estimate based on Eq. 13 for typical experimental parameters in sodium (d — 10 fim, 
peak density n = 1.5 x 10 14 cm -3 , a = 2.75 nm) yields v c ~ 1.0 mm/s, close to 
the experimental observations. However, Eq. 13 depends only weakly on the speed 
of sound, through the logarithmic dependence on the healing length £. In contrast, 
our measurements show an approximate proportionality to the sound velocity [10], 
suggesting that vortex nucleation, determines the onset of dissipation. 

Time-dependent simulations of the Gross-Pitaevskii equation show the formation 
of vortex line pairs, above a critical velocity which is close to the observed value 
[33] . Several authors have emphasized the role of locally supersonic flow around the 
laser beam in the nucleation of vortices. [33-35]. In one theoretical model [34], the 
vortices are emitted periodically at a rate that increases with velocity, and reduce 
the pressure gradient across the object. The predicted heating rate [34,36] is in 
rough agreement with the data. Moreover, this model also predicts that the slope 



of the asymmetry should increase at lower density, in accord with our observations. 



V AMPLIFICATION OF LIGHT AND ATOMS 

Spontaneous light scattering can be stimulated when the atomic recoil state is 
already populated (the N q term in Eq. 6). We have explored this process in our 
studies of superradiance [12], phase-coherent atom amplification [13] and optical 
amplification [14]. 

In all these experiments, the condensate was illuminated with a laser beam (mode 
k, also called "dressing beam"). A condensate atom scatters a photon from the 
laser beam into another mode and receives the corresponding recoil momentum and 
energy. Injection of atoms or light turns this spontaneous process into a stimulated 
process and realizes an amplifier for either atoms or light. If atoms are injected, 
they form a matter wave grating (an interference pattern with the condensate at 
rest) and this grating diffracts light. The diffraction transfers recoil momentum 
and energy to the atoms, which results in a growth of the grating and therefore the 
number of atoms in the recoil mode — this is the intuitive picture for atom gain. 

Eq. 6 describes gain for atoms. In the limit of an empty mode for the scattered 
light {n k _ q = 0), one obtains 

2n 

W + = —\C\ 2 N n k (N q + l)5(E k - E k _ q - hoof) (17) 

For the high momentum transfers considered here (on the order of the photon 
recoil momentum), S(q) = 1. Each scattering event which transfers momentum 
q to the condensate, generates a recoiling atom — therefore the scattering rate W + 
integrated over all final states for the scattered photon gives the growth rate N q 
for the recoiling atoms: 

N q = {G q -Y 2A )N q (18) 



with the gain coefficient 



G - RN ^An (19) 
q 8tt/3 q y ' 

Here R is the rate for single-atom Rayleigh scattering which is proportional to 
the pump light intensity, No the number of atoms in the condensate at rest, 9 q 
the angle between the polarization of the incident light and the direction of the 
scattered light, and fl q the phase-matching solid angle for scattering into mode q. 
In addition, a loss term T 2 , q was included which describes the decoherence rate of 
the matter-wave grating and determines the threshold for exponential growth. It 
represents the linewidth of the two-photon process generating recoil atoms in mode 
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FIGURE 9. Input-output characteristic of the matter-wave amplifier. (a-c) Typical 
time-of-flight absorption images demonstrating matter wave amplification. The output of the 
seeded amplifier (c) is clearly visible, whereas no recoiling atoms are discernible in the case with- 
out amplification (a) or amplification without the input (b). The size of the images is 2.8 mm x 
2.3 mm. (d) Output of the amplifier as a function of the number of atoms at the input. A straight 
line fit shows a number gain of 30. Reprinted by permission from Nature, Ref. [13], copyright 
1999 Macmillan Magazines Ltd. 



FIGURE 10. Experimental scheme for observing phase coherent matter wave amplification. A 
small-amplitude matter wave was split off the condensate by applying a pulse of two off-resonant 
laser beams (Bragg pulse). This input matter wave was amplified by passing it through the 
condensate pumped by a laser beam. The coherence of the amplified wave was verified by ob- 
serving its interference with a reference matter wave, which was produced by applying a second 
(reference) Bragg pulse to the condensate. The interference signal was observed after 35 ms of 
ballistic expansion. The fringes on the right side show the interference between the amplified 
input and the reference matter wave. Reprinted by permission from Nature, Ref. [13], copyright 
1999 Macmillan Magazines Ltd. 
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Fig. 9 shows the input-output characteristics of the amplifier. The number of 
input atoms was below the detection limit of our absorption imaging (Fig. 9a). The 
amplification pulse alone, although above the threshold for superradiance [12], did 
not generate a discernible signal of atoms in the recoil mode (Fig. 9b). When the 
weak input matter wave was added, the amplified signal was clearly visible (Fig. 9c). 
The gain was controlled by the intensity of the pump pulse (see Eq. 19) and typically 
varied between 10 and 100. Fig. 9d shows the observed linear relationship between 
the atom numbers in the input and the amplified output with a number gain of 30. 

The phase of the amplified matter wave was determined with an interferometric 
technique. For this, a reference matter wave was split off the condensate in the 
same way as the first (input) wave (see Fig. 10). The phase of the reference matter 
wave was scanned by shifting the phase of the radio-frequency signal that drove the 
acousto-optic modulator generating the axial Bragg beam. We then observed the 
interference between the reference and the amplified matter waves by measuring 
the number of atoms in the recoil mode. 

The atom amplification is described by the Hamiltonian (Eq. 1) as a four-wave 
mixing process between two electromagnetic fields and two Schrodinger fields. The 
symmetry between light and atoms indicates that a dressed condensate should 
not only amplify injected atoms, but also injected light. In Ref. [16] we have 
discussed that matter wave gain and optical gain emerge as two limiting cases, 
depending on whether the atomic population N q or occupation of the optical mode 
rik-q dominates the bosonic stimulation term in Eq. 7. In our experiments, the 
recoiling atoms move ten orders of magnitude more slowly than the scattered light. 
Therefore, we always have N q ^> n^q and n k ^ q < 1 and optical stimulation should 
play no role. Still, as we want to discuss now, the dressed condensate can act as 
an amplifier for light. 

The physical picture behind the optical gain is as follows: If a very weak probe 
beam is injected into the dressed condensate, it acts together with the dressing 
beam as a pair of Bragg beams and creates recoiling atoms. Those recoiling atoms 
move out of the condensate (or decohere) on a time scale T^" 1 which is the inverse 
of the linewidth of the Bragg transition. In steady state, the number of recoiling 
atoms N q in the volume of the condensate is proportional to the intensity of the 
probe light. Those recoiling atoms interfere with the condensate at rest and form 
a diffraction grating which diffracts the dressing beam into the path of the probe 
light resulting in amplification of the probe light (Fig. 11). 

An expression of the gain can be derived in analogy to a fully inverted two- 
level system with dipole coupling which would have a gain cross section of 67rA 2 
for radiation with wavelength A(= 2nX). For the Raman-type system in Fig. lib, 
the gain is reduced by the excited state fraction, R/Y\ (where R is the Rayleigh 
scattering rate for the dressing beam and Ti is the linewidth of the single-photon 
atomic resonance) and increased by Ti/T 2 , the ratio of the linewidths of the single- 
photon and two-photon Bragg resonances. Thus the expected cross-section for gain 
is 
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FIGURE 11. Amplification of light and atoms by off- resonant light scattering, (a) The fun- 
damental process is the absorption of a photon from the "dressing" beam by an atom in the 
condensate (state |1)), which is transferred to a recoil state (state |2)) by emitting a photon into 
the probe field. The intensity in the probe light field was monitored by a photomultiplicr. (b) The 
two-photon Raman-type transition between two motional states (|1), |2)) gives rise to a narrow 
resonance, (c) The dressed condensate is the upper state (|1')) of a two-level system, and decays 
to the lower state (recoil state of atoms, |2)) by emitting a photon. Figure is taken from Ref. [14] 



<Tgain = 67rA 2 — • (20) 
1 2 

We observed quasi-steady state gain by a factor of up to three. This optical 
gain has a narrow bandwidth due to the long coherence time of a condensate. The 
gain represents the imaginary part of the complex index of refraction. A sharp 
peak in the gain implies a steep dispersive shape for the real part of the index of 
refraction n(uj). This resulted in an extremely slow group velocity for the amplified 
light. Fig. 12 shows that light pulses were delayed by about 20 /is across the 20 /zm 
wide condensate corresponding to a group velocity of 1 m/s. This is one order of 
magnitude slower than any value reported previously (see Ref. [37] and references 
therein). 

The (amplitude) gain g for the probe light is related to the matter wave gain G 

g=l + ^l = 1+ C (21) 
^ J- 2 

where n is the condensate density and / its length. However, when the gain G 
is above the threshold for superradiance, G > T 2 (Eq. 18) the optical gain should 
diverge: a single recoiling atom created by the probe light and dressing light is 
exponentially amplified and creates a huge matter wave grating which will diffract 
the dressing light into the probe light path, thus amplifying the probe light by a 
large factor. In order to describe this non-linear feedback, one has to use coupled 
equations for optical gain and atom gain [14]. As a result, one has to replace the 
loss rate T 2 in Eq. 21 by the dynamic rate T 2 — G (see Eq. 18) 



which agrees with Eq. (21) in the low-intensity limit. By raising the gain over 
the threshold we could map out the transition from single-atom gain to collective 
gain [14]. The expansion g = 1 + (G/r 2 ) + (G/T 2 ) 2 + . . . shows the transition 
from (linear) single-atom gain to (non-linear) collective gain and illustrates that 
the dressed condensate is a clean model system for discussing linear and non-linear 
behavior, optical and atom-optical properties and their interplay. 
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FIGURE 12. Pulse delay due to light amplification, (a) About 20 /is delay was observed when 
a Gaussian pulse of about 140 [is width and 0.11 mW/cm 2 peak intensity was sent through the 
dressed condensate (bottom trace). The top trace is a reference taken without the dressed con- 
densate. Solid curves are Gaussian fits to guide the eyes, (b) The observed delay was proportional 
to the logarithm of the observed gain. Figure is taken from Ref. [14] 



VI DISCUSSION 

This paper has summarized our recent experiments on Bose-Einstein conden- 
sation with the unifying theme of suppression and enhancement. Suppression of 
scattering or dissipation can arise for two different reasons. The phonon and vortex 
nature of the collective excitations together with energy and momentum conserva- 
tion allow dissipation only above a critical velocity. In addition, one has to consider 
the dynamics of the excitation process. For microscopic particles, this is reflected in 
the matrix element S(q) which characterizes how easily can the condensate absorb 
momentum. For macroscopic motion, it is reflected in a critical velocity for vortex 
nucleation. Scattering processes are also enhanced by the population in the final 
states (bosonic stimulation). Optical stimulation was used in Bragg scattering, and 
matter wave stimulation led to superradiance and matter wave amplification. 

In closing, let us note that there are some subtleties which go beyond the sim- 
ple picture using rate equations and occupation numbers. A condensate in its 
ground state is in a coherent superposition state of the zero-momentum state with 
correlated pairs with momenta ±q (the quantum depletion). The population in 
the quantum depletion can cause bosonic stimulation of spontaneous emission [38]. 
However, for a scattering situation, there are two bosonically enhanced pathways 
which destructively interfere (causing S(q) < 1, Sect. II). Therefore, the concept 
of bosonic stimulation can be applied to the quantum depletion, but with caution 
[38]. 



Finally, the matter wave amplification described in Sect. V can be regarded as 
the self-amplification of a density modulation (caused by matter wave interference 
between the condensate and the amplified recoiling atoms). A similar amplification 
can happen for a density modulation in a fermionic gas [39], but, the coherence time 
of such a density modulation is generally much shorter than in a condensate. This 
example shows that a more general description of stimulation and amplification has 
to address the symmetry and coherence of the prepared state. Bosonic quantum- 
degeneracy is sufficient, but not necessary for stimulated scattering [39]. 

We are grateful to Axel Gorlitz for helpful discussions. This work was supported 
by NSF, ONR, ARO, NASA, and the David and Lucile Packard Foundation. 
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